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ABSTRACT 

SN 2006jc is a peculiar supernova (SN), in which the formation of dust has 
been confirmed at an early epoch of ~50 days after the explosion. We investigate 
the possibility of such an earlier formation of dust grains in the expanding ejecta 
of SN 2006jc, applying the Type lb SN model that is developed to reproduce the 
observed light curve. We find that the rapid decrease of the gas temperature in 
SN 2006jc enables the condensation of C grains in the C-rich layer at 40-60 days 
after the explosion, which is followed by the condensation of silicate and oxide 
grains until ~200 days. The average radius of each grain species is confined to 
be less than 0.01 /xm due to the low gas density at the condensation time. The 
calculated total dust mass reaches ~1.5 M Q , of which C dust shares 0.7 M & . 
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On the other hand, based on the calculated dust temperature, we show that the 
dust species and mass evaluated to reproduce the spectral energy distribution 
observed by AKARI and MAGNUM at day 200 are different from those obtained 
by the dust formation calculations; the dust species contributing to the observed 
flux are hot C and FeS grains with masses of 5.6 x 10 -4 M & and 2.0 x 10~ 3 M & , 
respectively, though we cannot defy the presence of a large amount of cold dust 
such as silicate and oxide grains up to 0.5 M . One of the physical processes 
responsible for the difference between calculated and evaluated masses of C and 
FeS grains could be considered to be the destruction of small-sized clusters by 
energetic photons and electrons prevailing within the ejecta at the earlier epoch. 

Subject headings: dust, extinction — supernovae: general — supernovae: indi- 
vidual (SN 2006jc) 

1. Introduction 

Supernova (SN) 2006jc is a peculiar Type lb SN (SN lb) in which the formation of 
dust has been observationally confirmed at a very early time of ~50 days after the peak 
brightness (Smith et al. 2008), which is more than a few hundred days earlier than those 
for the dust-forming SNe observed so far. Evidence for dust formation in SN 2006jc comes 
from an increase of the red to near-infrared (NIR) continuum (Di Carlo et al. 2007) and 
simultaneous emergence of blueshifted narrow He I emission lines (Smith et al. 2008). Smith 
et al. (2008) have proposed that the detected dust does not condense in the ejecta but forms 
in the dense postshock gas swept up by the forward and reverse shocks. Although X-ray 
observations of SWIFT and Chandra (Immler et al. 2008) suggested the presence of the 
dense circumstellar (CS) shell ejected by the outburst similar to those seen in luminous blue 
variables (LBVs) two years prior to the explosion (Nakano et al. 2006; Pastorello et al. 
2007), the postshock CS gas density estimated from the X-ray light curve (Tominaga et al. 
2007) is not high enough for dust grains to condense. Thus, it is reasonable to consider that 
the appearance of dust toward SN 2006jc should be the outcome of ongoing dust formation 
in the expanding ejecta of the SN. 

Several pieces of evidence for dust formation in the SN ejecta have been reported for 
Type II SNe (SNe II): SN 1987A (Lucy et al. 1989; Whitelock et al. 1989; Meikle et al. 
1993; Wooden et al. 1993; Colgan et al. 1994), SN 1998S (Gerardy et al. 2000; Pozzo 
et al. 2004), SN 1999em (Elmhamdi et al. 2003), and SN 2003gd (Sugerman et al. 2006; 
Meikle et al. 2007). For these SNe II except for SN 1998S, the onset of dust formation is 
estimated to be later than 400 days after explosions. Type Iln SN 1998S with the relatively 
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low-mass hydrogen envelope (e.g., Liu et al. 2000; Fransson et al. 2005) exhibits signatures 
of dust condensation around 230 days (Gerardy 2000). In addition, SN 19901 is the first 
SN lb signifying the ongoing dust formation, in which dust formation is observed at ~ 230 
days (Elmhamdi et al. 2004). On the other hand, the signature of dust formation in SN 
Ic has not been recorded so far. Although theoretical studies have shown feasibility of dust 
formation in SNe II (Kozasa et al. 1989, 1991; Todini & Ferrara 2001; Nozawa et al. 2003), 
the possibility of dust formation in SNe Ib/c has never been explored to date. Hence, it is 
an important subject to be pursued how the formation process of dust in the ejecta depends 
on the type of SNe. 

Since the progenitor stars of SNe Ib/c have lost most of the hydrogen/helium envelopes 
before the explosion, their ejected masses are smaller, and the expansion velocities are sig- 
nificantly higher than SNe II. This leads to a lower density of gas in the ejecta, and the gas 
temperature drops down more quickly than those in typical SNe II. This allows us to expect 
earlier dust formation in SNe Ib/c than SNe II. On the other hand, the lower gas density 
may even result in unsuccessful dust formation. Furthermore, in the early epoch, energetic 
photons and electrons generated by the Compton degradation of 7-rays from the decay of 
radioactive elements such as 56 Ni and 56 Co prevail throughout the ejecta, and may affect 
the formation process of dust grains. SN 2006jc is an ideal laboratory for examining the 
process of dust formation in SNe lb because a copious amount of observational data enable 
us to compare with theoretical models. In this paper, we investigate, for the first time, the 
possibility of dust formation in the expanding ejecta of SN lb, based on the SN model that 
can well reproduce the light curve of SN 2006jc (Tominaga et al. 2007). 

The paper is organized as follows. In § 2 we summarize the observational evidence for 
dust formation in SN 2006jc and its interpretation, which are compared with the calculations 
of dust formation in the later sections. In § 3 we investigate the dust formation in the ejecta 
and show that it is possible for dust grains to condense in the ejecta of SN 2006jc at very 
early times of about 50 days after explosion. In § 4 we calculate the temperature of possible 
condensates obtained in § 3 and evaluate the amount of dust grains contributing to the 
spectral energy distribution observed with AKARI and MAGNUM at day 200. Formation 
processes of dust grains in SN 2006jc are discussed in § 5, focusing on the interpretation 
of the difference between dust mass obtained by the dust formation calculation and that 
evaluated from the spectral fitting. The summary is presented in § 6. In this paper, we 
assume that the explosion of SN 2006jc has occurred at 15 days before its discovery. 
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2. Observations of Dust Formation in SN 2006jc and Its Interpretation 

2.1. Evidence for dust formation 

The first verification of the presence of dust in SN 2006jc comes from the re-brightening 
at J, H, and K bands around day 50 (Arkharov et al. 2006; Di Carlo et al. 2007). In 
addition, Di Carlo et al. (2007) have noted a concomitant steep drop of the light curve in 
the optical bands starting from 70 days, suggestive of strong extinction by the newly formed 
dust. It is also shown that this rapid decline of the optical light curve continues up to ~120 
days (Tominaga et al. 2007; Anupama et al. 2008; Kawabata et al. 2008; Mattila et al. 
2008). 

Definitive evidence for dust formation has been reported by Smith et al. (2008); the 
asymmetry of narrow He I emission lines increased concurrently with increasing the red 
to NIR continuum between 65 and 120 days. The increasing red/NIR continuum is inter- 
preted as the thermal radiation from the newly formed carbonaceous dust with temperature 
of ~1600 K, considering its high temperature. At day 140, the excess of this continuum 
disappeared, which indicates that the dust grains were completely destroyed or the dust 
temperature became so low that its thermal emission was shifted to longer wavelengths. 
The increasing asymmetry of narrow He I emission lines tracing the shocked CS gas was 
caused by the dust obscuration of the redshifted side, whereas broad emission lines trac- 
ing the SN ejecta such as Ca II and O I faded entirely. The blueshifts of He I emission 
lines persisted up to day 180 (Mattila et al. 2008), asserting that the newly formed carbon 
grains were not destroyed but cooled down, which is supported by the NIR observations 
with AKARI (Sakon et al. 2007) and MAGNUM (Minezaki et al. 2007) showing that the 
temperature of the carbon grains dropped down to ~800 K at day 215. 

The time evolution of the optical light curve and the He I emission line profile during 
50-140 days strongly suggests the ongoing formation of dust in SN 2006jc and reliably rules 
out as origins of the red/NIR continuum the scenarios involving pre-existing dust in the CS 
medium (CSM), such as the IR light-echo (Smith et al. 2008, Mattila et al. 2008) and the 
heating of the CS dust through the collision with the forward shock. However, as suggested 
by Sakon et al. (2007) and Mattila et al (2008), the mid-IR spectrum observed at later epoch 
(> 200 days) can be attributed to the low-temperature CS dust, which is discussed in § 4.2. 
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2.2. The site of dust formation 

Smith et al. (2008) have supposed that dust grains are formed in the dense shell swept up 
by shocks but not in the freely expanding ejecta, based on (1) the geometrical considerations 
for explaining the blueshifts of narrow He I lines and the overall fading of broad Ca II and 
O I lines and (2) the time evolution of dust temperature that remains nearly constant for 
about a month and decreases after then; they considered that dust can continuously form at 
temperature of ~1600 K in the rapidly cooling regions behind the forward shock rather than 
the reverse shock, in the course of the passage of the forward shock through the dense CS 
shell created by the LBV-like eruption two years prior to the SN explosion. This scenario is 
supported by a flux ratio of He I A7065/A5876 close to unity and the He II A4686 emission 
feature appeared during the putative time of dust formation. The detection of these He 
lines are expected to indicate the presence of the postshock gas of ~3 x 10 -14 g cm -3 that 
is comparable to the critical density for carbon grain condensation (e.g., Clayton 1979), 
although Mattila et al. (2008) can not find the sign of the He II line in their spectrum 
during 88-132 days. 

Here we investigate whether dust formation in the postshock gas around SN 2006jc is 
possible or not, referring to the X-ray observations and the hydrodynamical calculations for 
the ejecta-CSM interaction. Observations with SWIFT and Chandra satellites (Immler et 
al. 2008) have shown that the X-ray luminosity is in the range of 10 39 -10 40 ergs s -1 from 
day 36 through day 180 with a gentle rise of the X-ray emission until 110 days. Immler 
et al. (2008) have estimated the density of the postshock gas shell to be ~4 x 10~ 17 g 
cm -3 , assuming that all of the measured X-rays at peak are emitted from the dense shocked 
shell. However, the time evolution of the X-ray luminosity indicates that the medium around 
SN 2006jc has a quite smooth density profile. Tominaga et al. (2007) have performed the 
hydrodynamic calculations to reproduce the observed X-ray light curve and found that the 
density profile of the preshock CSM should be po = 2.75 x 10~ 19 g cm -3 for r < 2.2 x 10 16 
cm and p = po(r/2.2 x 10 16 cm)~ 6 g cm" 3 for r > 2.2 x 10 16 cm. Their result has shown that 
the gas density swept up by the forward and reverse shocks is < 10~ 17 g cm" 3 at 50 days 
(see Figure 10 in Tominaga et al. 2007), which is too low for dust grains to nucleate. Thus, 
on the basis of the hydrodynamical simulations, we conclude that dust grains are not likely 
to be formed in the postshock CS gas. 

On the other hand, it should be addressed that there are some studies that advocate 
formation of dust in the dense CS shell produced by the LBV-like outburst before two years 
(Di Carlo et al. 2007; Mattila et al. 2008). Mattila et al. (2008) have shown that the 
observed X-ray luminosity is produced by the CS shell with the preshock density of (2.5-4) 
x 10 -15 g cm -3 in aid of their calculations of shock evolution and that the NIR light curve 
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is reproduced by the IR light echo of the dust newly formed in the dense, cooling CS shell. 
However, formation process of dust heavily depends on the cooling history of the gas, and 
whether dust formation can be realized in the rapidly cooling postshock gas should be fully 
explored in the future work. In the following sections, we shall persue the possibility of dust 
formation in the expanding ejecta of SN 2006jc. 

3. Dust Formation in the Ejecta of SN 2006jc 

3.1. Model of SN 2006jc and calculation of dust formation 

The time at which dust condenses in the SN ejecta and the size and number density 
of the newly formed dust are sensitive to the time evolution of the gas temperature and 
density (Kozasa et al. 1989; Todini & Ferrara 2001; Nozawa et al. 2003), given the elemental 
composition. In the calculation of dust formation, we apply the results of hydrodynamics and 
nucleosynthesis calculations for the SN lb model well reproducing the observed bolometric 
light curve of SN 2006jc (Tominaga et al. 2007); the ejecta mass is 4.9 M & with the explosion 
energy of i? ex /10 51 ergs = E 5 i = 10 and the ejected 56 Ni mass of 0.22 M , where 56 Ni is 
uniformly distributed in mass within the He core. The pre-explosion mass of the progenitor 
is 6.9 Mq. The elemental composition in the ejecta has an original onion-like composition, 
except for 56 Ni. The time evolution of the gas temperature in the freely expanding ejecta 
was calculated by solving the radiative transfer and the energy equations, taking account of 
the energy deposition from 56 Ni and 56 Co. 

Figure 1 shows the structures of the gas temperature (Fig. la) and density (Fig. 16) 
within the He core of the SN 2006jc model at 50 days (red solid lines) and 200 days (red 
dashed lines) after the explosion. For comparison, shown are those for the SN II model with 
the progenitor mass M pr = 25 M & and E§\ = 1 (blue solid and dashed lines at 50 and 200 
days, respectively; Umeda & Nomoto 2002), which has a He core of 5.5 M Q comparable to 
that of SN 2006jc. The gas temperature in SN 2006jc decreases much earlier, and the density 
at a given time is more than three orders of magnitude lower than that in the SN II. This is 
because the expansion velocity in the He core of SN 2006jc is considerably high, reflecting 
the absence of the hydrogen envelope and a ten times higher explosion energy than SNe II 
with E51 = 1. The gas temperature in SN 2006jc drops down to ~2000 K at 50 days and 
~1000 K at 200 days. Thus, we can expect early formation of dust grains between 50 days 
and 200 days, since a typical dust condensation temperature ranges from 1000 K to 2000 K. 

The dust formation calculation is performed by applying a theory of non-steady state 
nucleation and grain growth described in Nozawa et al. (2003). In the calculations, we 
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assume the sticking probability of a s = 1, and that the temperature of a small cluster 
consisting of up to a few tens atoms is the same as the gas temperature^ 

3.2. Results of dust formation calculations 

Figures 2 a and 2b show the results of calculation. Figure 2a indicates the condensation 
time of the newly formed dust defined as the time when the nucleation rate reaches the 
maximum (see Nozawa et al. 2003). It can be seen that in the outermost C-rich layer of 
the ejecta (see Fig. 3c in Tominaga et al. 2007), C grains condense at very early times of 
40-60 days. This is due to the rapid decrease of the gas temperature in the ejecta as well as 
the relatively high condensation temperature of C grains. It should be emphasized that the 
condensation times of C grains are in good agreement with the onset time of dust formation 
observed in SN 2006jc. Furthermore, the condensation of C grains can explain both the 
blueshift of He I lines and the entire fading of Ca II and O I lines because the formation site 
is in the region between the He-rich CSM and the O-rich layer in the ejecta, though they 
cannot absorb the zero-velocity component of narrow He I emission lines efficiently. 

After C grains are formed, AI2O3, Mg 2 Si04, MgSiOs, SiC>2, Fe3C>4, MgO, and FeO grains 
form in the O-rich layer at 90-150 days in this order, and FeS and Si grains condense in the 
Si-S-rich layer at ~200 days. Since 56 Ni is uniformly mixed throughout the ejecta, Fe3C>4 and 
FeO grains are produced instead of Fe grains; Fe grains are expected to form in the innermost 
Fe-Ni core unless 56 Ni intrudes into the outer layers (Nozawa et al. 2003). However, the 
mixing of elements is likely to be at the knotty level rather than at the microscopic atomic 
level, as suggested by the observations of the Cas A supernova remnant (Douvion et al. 2001; 
Ennis et al. 2006; Rho et al. 2008). If 56 Ni is not mixed at the atomic level, the formation 
of Fe grains is expected as well. Nevertheless, Fe grains, which have a high energy barrier 
for nucleation and a resulting low condensation temperature (<800K), cannot form in the 
ejecta considered here because the gas density is too low at the time when the gas cools down 
to 800K (~250 days). Also, Fe atoms originally existing in the O-rich layer are not locked 
into grains significantly (condensation efficiency of <10~ 7 ) owing to their small abundances. 
Considering these uncertainties, we rule out Fes04 and FeO grains as the newly formed dust 
species. 

As is expected from the low gas density as well as the rapid cooling of gas in the ejecta, 
the average radius of each dust species given in Figure 2b is considerably small, being less 



1 There is a typographical error in the equation (3) for the nucleation rate in Nozawa et al. (2003); 
replace cij with Cy. 
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than 0.01 /im for all grain species. On the other hand, the condensation efficiency of dust 
is almost unity, and the amount of dust formed is relatively large. The total mass of newly 
formed dust is 1.45 M , and the mass of C grains (0.7 M Q ) accounts for about a half of the 
total mass. The mass of each dust species M±j obtained by the dust formation calculation 
is summarized in Table 1. Other major grain species in mass are Si02 (0.23 M ), Si (0.2 
M ), and MgSi0 3 (0.16 M ). 

It should be noted here that the dust formation process as demonstrated above is not 
influenced by the interaction between the ejecta and the CSM. Our results show that dust 
grains can not condense in the region of 0.005 M from the outer edge of the ejecta where 
the gas density is too low. Although the ejecta-CSM interaction causes the reverse shock 
penetrating into the ejecta, the hydrodynamic simulation shows that the reverse shock sweeps 
up only ~0.011 M of the outermost region in the ejecta by 50 days and ~0.053 M & by 200 
days (see Fig. 10 in Tominaga et al. 2007). Thus, the reverse shock sweeps up 0.006 M 
from the outer edge of the dust formation region at day 50, and 0.048 M & at day 200. The 
mass of C grains within the region swept up by 200 days is only 1.3 x 10~ 3 M , which is 
much smaller than 0.7 M at the time of formation, thus almost all of C grains can form 
in the freely expanding ejecta. Since the radii of C grains formed in the outermost region 
are smaller than 10 A (see Fig. 2b), the timescale of dust destruction by sputtering is 
estimated to be <6 hours (Nozawa et al. 2006) in the shocked gas with density of ~10~ 18 
g cm -3 (Tominaga et al. 2007). Hence, thermal emission from shock-heated grains can be 
negligible. 

4. Evaluation of Dust Temperature and Mass in the Ejecta 

The NIR to mid-IR (MIR) observations with AKARI (Sakon et al. 2007) and MAGNUM 
(Minezaki et al. 2007) have detected thermal emission from dust in SN 2006jc at ~200 days. 
Sakon et al. (2007) have concluded that the obtained NIR spectrum is fitted by amorphous 
carbon grains of 800 K with a mass of 6.9 x 10~ 5 M , assuming optically thin thermal 
radiation. They have also shown that the MIR excess emission over the amorphous carbon 
grains of 800 K can be explained either by silicate and silica grains of 700-800 K with ~10~ 4 
Mq or amorphous carbon grains of 320 K with 2.7 x 10~ 3 M . In this section, referring to 
the grains species calculated in § 3 as possible condensates in the ejecta and deriving the 
dust temperature, we evaluate the dust mass by reproducing the spectral energy distribution 
observed at day 200 with AKARI and MAGNUM. 
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4.1. Temperature of dust 

The equilibrium temperature of dust in the ejecta is determined by the balance between 
heating and cooling through radiation and collision with gas. In order to calculate the heating 
due to the absorption of ambient radiation, we evaluate the flux F\(r,t) at a time t and a 
position r within the ejecta, assuming that the total flux F(r,t) calculated by Tominaga et 
al (2007) radiates as a blackbody with temperature Tbb(?", t). 

Given the gas temperature T gas (r, t) and number density n gas (r, t) in the ejecta, the 
equilibrium temperature of dust T d (r, t) at a given time and at a given position is determined 
by implicitly solving the equation 

4vraV B T d (r,t) 4 (g A (a,T d )) = f° ' ira 2 Q x (a)B x (T BB )d\ 

vb+bb Jo 

\27r/im H / \27r/im H / 



2 %as ( c v + -k 



T 



where a is the radius of dust, Ob is the Stefan-Boltzmann constant, (Q\(a, T)) is the Planck- 
mean of absorption coefficient, B\(T) is the Plank function, k is the Boltzmann constant, 
c v is the heat capacity of a gas molecule (c v = 3k/ 2 for a monatomic molecule), \i is the 
mean molecular weight, and mn is the atomic mass unit. Here we consider that the gas 
colliding with dust is completely accommodated and is ejected from the surface with the 
same temperature as the dust (see Landau & Pitaevski 1981). In the calculations we assume 
^bb — 5000 K throughout the ejecta, regardless of time, and a = 0.01 /mi. Since in the 
ejecta considered here, the total flux at day 200 is F ~ 10 6 erg s _1 cm -2 almost uniformly 
over the He core, the energy transfered through radiative process is much larger than that 
by collision with gas, and thus we can neglect the collisional heating and cooling. Note that 
in this case the calculated dust temperature is independent of dust size because Q\(a)/a 
does not depend on the radius for a < 0.05 /mi. The references for the optical constants of 
each dust species used in the calculations are given in Table 1. 

We also note that the effect of stochastic heating can be negligible as long as we assume 
a blackbody radiation with Tbb ~ 5000 K under the given total flux. We have confirm that 
for C grains of 10 (20) A, stochastic heating only makes the dust temperature distribute sym- 
metrically around the equilibrium temperature with the temperature width of ~±100 (30) 
K. Accordingly, the resulting emissivity of 10 A-sized C grains does not change significantly 
at wavelengths longer than 2 /mi, compared with that derived for equilibrium temperature; 
its difference is 10 % at 2 /tm and is less than 2.5 % at > 3 /mi. Therefore, stochastic heating 
does not affect the results of the IR spectrum energy distribution calculated in § 4.2. 



Figure 3 shows the temperatures of all dust species at 200 days (solid lines), along with 
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the temperature of C grains at 60 days (dashed lines) a The temperature of the C grains 
is 1200-1700 K at 60 days and is higher for those formed in the inner region. Note that 
this high temperature is in good agreement with the temperature of dust responsible for 
the red/NIR continuum appearing in SN 2006jc two months after the explosion (Smith et 
al. 2008). At 200 days, the C grains cool down to 500-600 K, which is also considered to 
be correspondent to the temperature of 800 K for amorphous carbon grains suggested by 
the observations with AKARI and MAGNUM (Sakon et al. 2007), although the calculated 
temperature is relatively low. Thus, the calculated temperature of C grains can reasonably 
explain the time evolution of the dust temperature indicated by the observations toward SN 
2006jc from day 50 through day 200. 

On the other hand, silicate grains have a temperature of only 100-200 K at day 200, 
reflecting much lower absorption efficiency at optical wavelengths than at IR wavelengths. 
This low temperature cannot explain the MIR excess emission over the amorphous carbon 
grains with 800 K observed by AKARI (Sakon et al. 2007), and we can conclude that neither 
silicate nor silica grains newly formed in the ejecta can be the carriers of the MIR excess 
emission. In addition, our result shows that the temperature of C grains spans the range from 
500 to 600 K. Therefore, if the MIR excess emission would be attributed to the amorphous 
carbon grains with 320 K, they might not be the grains formed in the ejecta but could be 
the pre-existing CS grains heated by the SN outburst, as suggested by Sakon et al. (2007). 

In the calculations we assumed the ambient radiation field as a blackbody with Tbb = 
5000 K. The calculated dust temperatures at day 200 is not sensitive to Tbb as long as Tbb 
is in the range of 5000 K to 6000 K; for Tbb = 6000 K, the temperature of MgSiC>3 and 
Mg2SiC>4 grains increases by ~50 K compared with that for Tbb = 5000 K, while the increase 
of temperature of the other grains species is at most 30 K. Here, it should be kept in mind 
that, being different from normal SNe lb, SN 2006jc shows the strange optical spectrum 
in the early times until day 120, with a bright blue continuum whose origin is not clear 
(Foley et al. 2007, Pastorello et al. 2007). This indicates that the spectrum of radiation at 
the early epoch within the ejecta of SN 2006jc is not also expected to be approximated by 
a blackbody radiation. Nevertheless, as shown above, we can reasonably explain the time 
evolution of temperature of C grains indicated by the observations. Thus, we adopt the 
results for T B b = 5000 in the following discussions. 



2 Note that in the calculation we exclude Si grains; the temperature of Si grains at the condensation time 
of ^200 days is about 300 K higher than the gas temperature. When we consider the non-LTE effect for 
dust formation, its condensation time is delayed to 240 days, as discussed in § 5.1, which is later than the 
observation date of AKARI and MAGNUM. 
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4.2. IR spectral energy distribution and dust mass 

In this subsection, we calculate thermal emission from dust grains based on the dust 
temperatures presented in Figure 3, and estimate the dust mass necessary for reproducing 
the spectral energy distribution (SED) at day 200 observed with AKARI and MAGNUM. 

Taking into account the self- absorption of dust, the observed flux density F£ hs of thermal 
emission from dust grains formed in the ejecta is evaluated by 

F? S = rE 47rr2m ^W^i?A(T d (r))exp[-(r A ( J R) - r x (r))]dr, (2) 

where D = 25.8 pc is the distance from observers (Pastorello et al. 2007), R is the outermost 
radius of the dust-forming region in the ejecta, m^j{r) is the mass of dust species j per volume 
at the position r, and K\,j = 3Qxj(a) / '4apj is the mass absorption coefficient, and the dust 
bulk density pj is taken from Nozawa et al. (2006). Note that k,\j does not depend on the 
size as long as a < 0.05 /im. The optical depth r A is given by 

rr 

T \( r ) = / y]™dj(r)K\jdr. (3) 



The magenta solid line in Figure 4 is the SED calculated by adopting the mass of each 
dust species M\j given in § 3.2, and the SED is dominated by thermal emission from C 
grains existing in the outermost ejecta as the result of the large optical depth due to their 
large mass. It is clear that the derived IR spectrum largely exceeds the observational data 
a t ^3 /im and cannot reproduce the observed SED. 

The NIR to MIR spectrum can be fitted using the derived grain temperatures because 
dust temperature does not depend on the mass of dust as can be seen from Equation (1), so 
the dust mass can be used as a free parameter in order to reproduce the observations. We 
calculate thermal emission from dust grains by taking the mass of each grain species at the 
position r as fjm<ij(r) and fit the SED to the observational data, where the value of a free 
parameter fj(< 1) is assumed to be constant for a given dust species. The best fitted result 
is displayed by the red solid line in Figure 4 and can well explain the observed NIR to MIR 
spectrum. The contribution of each grains species is depicted by the dashed lines, for which 
we adopt the masses of C and FeS grains and the upper mass limits of the other dust species 
given as M 2 j in Table 1. Here the mass of C grains is 5.6 x 10~ 4 M , reduced by a factor 
of 1250 in comparison with the mass obtained by the dust formation calculation. FeS grains 
can also contribute to the NIR spectrum because of its high temperature (~900 K), and its 
mass is reduced to 2 x 10~ 3 M Q (fj = 0.03). On the other hand, silicates and oxide grains 
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do not contribute to the NIR spectrum because their temperatures are too low (100-200 K). 
Furthermore, the large optical depth of MgSiC>3 and Mg 2 Si04 at MIR wavelengths prevents 
us from constraining the mass of these grain species. Thus, we can not defy the existence of 
a large amount of cold dust reaching up to ~0.5 M Q in the ejecta from the spectral fitting, 
and the dust mass derived under the optically thin assumption should be considered to be 
the lower limit (see also Meikle et al. 2007). 

We also investigate the effect of the presence of amorphous carbon grains in the CSM 
on the observed spectrum, which have been suggested to contribute to the MIR spectrum 
at >200 days (Sakon et al. 2007; Mattila et al. 2008). In the calculations, we take the 
temperature of CS carbon grains to be 320 K (Sakon et al. 2007). If their mass is smaller 
than 2.7 x 10 -4 M , the best-fitted spectrum in Figure 4 is not influenced by the pre-existing 
C grains. Further increase of CS dust mass can also reasonably reproduce the observed IR 
spectrum as seen from Figure 5, where the mass of newly formed and pre-existing C dust 
is 2.24 x 10~ 4 M Q and 1.4 x 10~ 3 M Q , respectively (the mass of other newly formed dust 
is the same as that in Figure 4). Therefore, the spectral fitting can allow the presence of 
low-temperature C grains in the CSM. 

The emission from CS grains can be the IR light-echo but not the thermal radiation 
from dust heated by the collision with the forward shock. This is readily demonstrated by 
considering the radius of dust cavity created by the evaporation of pre-existing CS grains 
due to strong radiation from a SN. Taking the evaporation temperature of C grains to be 
1800 K and using the SN peak luminosity of 10 42,7 ergs s -1 with the blackbody temperature 
of ~6000 K (Tominaga et al. 2007), we estimate the radius of dust cavity to be ~ 3 x 10 16 
cm. Because the estimated radius of dust cavity is comparable with the radius of the forward 
shock at ~200 days, the forward shock could have not swept a significant amount of CS dust. 
If we consider the radiation from the shock breakout as the heating source, which is expected 
to be brighter and have a higher blackbody temperature, the dust cavity may be extended 
so far that any CS dust can not be swept up by the forward shock at ~200 days. Note that 
in this spectral fitting, the mass of newly formed C grains depends on the temperature and 
mass of CS dust, and in the following discussion, we refer to the mass of C grains derived 
by the spectral fitting without CS dust in Figure 4. 
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5. Discussion 

5.1. Mass of dust formed in the ejecta 

As shown in § 4.2, the mass of C grains (0.7 M ) obtained by the dust formation 
calculation is three orders of magnitude larger than that (5.6 x 10 -4 M Q ) estimated from 
the spectral fitting to the observations. One of the reasons for this difference is considered 
to be the clumpy structure in the ejecta. In order to evaluate the dust mass necessary for 
reproducing the observed spectrum of SN 2003gd, Sugerman et al. (2006) performed the 
radiative transfer calculation taking account of the clumpy structure and found an increase 
in the estimated mass of dust in the ejecta by about one order of magnitude, compared with 
the mass derived for the uniform distribution of dust. Nevertheless, the consideration of dust 
clumping does not seem to be sufficient to overcome the large gap. Thus, the difference in 
dust mass should stem from the process of dust formation in the ejecta. 

One of the processes influencing on dust formation is the non-local thermal equilibrium 
(non-LTE) effect. As seen from Figures 1 and 3, the dust temperature is not the same as 
the gas temperature; at 200 days the temperature of C grains (500-600 K) and silicate and 
oxide grains (100-200 K) is lower than the gas temperature with 700-1000 K. In the ejecta 
with strong radiation field, the temperature of small clusters is expected to be higher than 
the gas temperature depending on the optical property. In this case, the non-LTE effect 
retards or hinders the formation of dust grains. Assuming that the optical properties of 
small carbon clusters are the same as those of small C grains, and applying a formula of 
nucleation rate taking into account the non-LTE effect by Kozasa et al (1996), we realized 
that the non-LTE effect does not significantly affect the formation of C grains; the non-LTE 
effect only retards the formation of C grains by less than ten days and does not affect the 
mass of C grains because the temperature of small carbon clusters is only 50 K higher than 
that of gas. Meanwhile, the non-LTE effect delays the formation of Si grains from ~200 days 
to ~240 days. This is the reason why we exclude Si grains in reproducing the observed SED 
at day 200 in § 4.2. 

Another process affecting the formation of dust in the ejecta is the energy deposition 
on small-sized clusters through the latent heat deposition at the condensation as well as the 
collisions with high energy photons and electrons. In the ejecta rich in condensible elements, 
the destruction of small-sized clusters by the deposition of latent heat could be significant 
unless the deposited energy is released by collisional and/or radiative processes efficiently 
Also, at an early epoch after the explosion, energetic photons and electrons generated from 
the decays of radioactive elements 56 Ni and 56 Co prevail abundantly throughout the ejecta. 
Since the nucleation rate is determined by the balance between growth and destruction rate 
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of small-sized clusters, the energy deposition makes small-sized clusters unstable against the 
growth and result reduce the nucleation rate. 

Although the nucleation rate can be reduced by increasing the critical cluster size as 
demonstrated by Bianchi & Schneider (2007), no reason is there to cut the number of seed 
nuclei by increasing the critical size arbitrarily, apart from the uncertainty inherent in the 
standard nucleation theory arising from the evaluation of chemical potential of small-sized 
clusters. Thus, we examine the effect of energy deposition for dust formation at early times 
observed in SN 2006jc by reducing the sticking probability a s for simplicity and comparing 
with the mass of dust evaluated in § 4.2. If we take a s ~ 3 x 10 -3 , we can obtain the 
mass of C grains of ~5 x 10 -4 M Q , while a s is ~0.3 to produce FeS grains of 2 x 10 -3 M Q . 
In this case, the condensation times of C and FeS grains is 40-67 days and 190-200 days, 
respectively. While the destruction rate of small-sized clusters depends on the properties of 
materials forming dust grains, the difference in the sticking probability may be caused by the 
collisions with energetic photons and electrons, reflecting the difference in the condensation 
time of C and FeS grains; energetic photons and electrons are more abundant at day 50 than 
at day 200. The details of these microscopic processes are beyond the scope of this paper 
and will be investigated in the future study on the dust formation in SNe. 



5.2. Early condensation of dust in SNe lb 

In § 3.2, we show that the early formation of dust around 50 days is possible in SN 
2006jc whose progenitor has lost most of He envelope (Foley et al. 2007; Tominaga et al. 
2007). However, other dust-forming SNe lb had not show such an early formation of dust; 
SN 19901 displayed the diagnostics of dust formation around 230 days (Elmhamdi et al. 

2004) , though the condensation time is relatively earlier than that of SN II. A peculiar Type 
lb SN 2005bf shows a rapid decline of the light curve from 50 days after the explosion, but 
did not exhibit the NIR brightening in contrast to the case of SN 2006jc (Tominaga et al. 

2005) . Maeda et al. (2007) have shown that the quickly faded light curve of SN 2005bf can 
be explained by the newly born magnetized neutron star at the center. On the other hand, 
they have demonstrated that the blueshifted profiles of O I and Ca II observed at ~270 days 
can be explained by the obscuration of dust in the ejecta. If the blueshift of emission lines 
in SN 2005bf is attributed to the newly formed dust grains, the time of dust formation is 
considered to be later than 100 days and not to be so early as in SN 2006jc, although the 
epoch of dust formation in SN 2005bf was not specified by the observations. 

The results of calculations revealed that the early formation of dust in SN 2006jc can 
be realized by the rapid decline of gas temperature, reflecting the large ratio of the explosion 
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energy (E 51 =10) to the ejected mass (M ej = 4.9 M ); the ratio E 51 /(M ci /M & ) ~ 2 is 
much greater than 0.125-0.25 for SN 2005bf (Tominaga et al. 2005) and ~0.3 for SN 19901 
(Elmhamdi et al. 2004). Unfortunately, no literature has reported the explosion energy of 
Type Iln SN 1998S, which had a time of dust condensation of ~230 days (Gerardy 2000). 
Using the ejecta mass estimated from the He core of 4-6 M (Fassia et al. 2001; Pooley et 
al. 2002) and hydrogen envelope of 0.1-2 M Q (Fassia et al. 2000; Anupama et al. 2001) 
and taking E 5i = 1 as a typical explosion energy, we can, however, deduce the ratio to be 
E 51 /(M ei /M Q ) = 0.125-0.25 for SN 1998S, which is similar to or slightly less than those for 
SN 2005bf and SN 19901. On the other hand, the ratio for the dust forming SNe II is of order 
0.1: 0.1-0.2 for SN 1987A (Shigeyama et al. 1988; Woosley 1988), 0.05-0.1 for SN 1999em 
(Elmhamdi et al. 2003), and ~0.12 for SN 2003gd (Hendry et al. 2005). Dust formation 
has not been yet observed on SNe Ic. However, the early dust formation in SNe Ic is also 
expected because SNe Ic may well have M e j similar to SNe lb. 

In conclusion, we can speculate that as the envelope of star at explosion is smaller and 
the explosion energy is higher, the condensation time of dust tends to be earlier. In order 
to prove these hypothesises and to investigate formation process of dust in SNe, we need 
further observational and theoretical studies on dust formation in various types of SNe. In 
particular, the UV to MIR observations at early times for SNe Ib/c are promising to clarify 
the formation process of dust because its condensation time is much earlier than that for 
SNe II, as is shown above. 

6. Summary 

We investigate formation of dust grains in the ejecta of SN 2006jc by applying the SN 
model that can nicely reproduce its light curve. We find that formation of dust in SN 2006jc 
is possible at early times, thanks to the rapid decrease of the gas temperature in the ejecta. 
The condensation of C grains is realized in the C-rich layer at very early times of ~50 days 
after the explosion. This is in good agreement with the time of the onset of dust formation 
indicated by the observations. In addition to C grains, the other grain species such as silicate 
and oxide grains are formed until ~200 day, with the average radii less than 0.01 /im. We 
suppose that formation of dust earlier than SNe II is common to the hydrogen-deficient SNe 
and that SN 2006jc is its extreme case. 

Furthermore, we calculate the dust temperature by considering the radiation field in the 
ejecta at the early epoch as a blackbody with Tbb = 5000 K. C grains have a temperature of 
~1600 K at 60 days and cools down to ~600 K at 200 days. This temperature of C grains 
can reasonably explain the time evolution of dust temperature suggested by the observations 
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toward SN 2006jc from day 50 through day 200. The temperatures of silicates and oxide 
grains are 100-200 K at day 200 and are too low to contribute to the MIR spectrum observed 
with AKARI. 

The comparison of the calculated SED with the observations of AKARI and MAGNUM 
shows that dust grains contributing to NIR/MIR flux are C and FeS grains whose masses 
are 5.6 x 10~ 4 M Q and 2 x 10~ 3 M , respectively. The masses of hot C and FeS grains 
obtained by the dust formation calculation are considerably larger than that estimated from 
the observations. This difference in dust mass may be explained by the destruction of small- 
sized clusters by the collisions with energetic photons and electrons. On the other hand, we 
can not constrain the mass of low-temperature silicate and oxide grains in the ejecta from 
the spectral fitting. The spectrum fitting also shows that the observed MIR spectrum can 
be attributed to pre-existing CS carbon grains with the mass up to ~10~ 3 M Q . In order 
to reveal the process of dust formation in the SN ejecta, it is essential to carry out UV to 
IR observations of various types of SNe. In particular, SNe Ib/c are the objects suitable 
for detecting signatures of dust formation by the optical to NIR observations because the 
condensation time of dust is much earlier than SNe II. 
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M r/ M He-core 

Fig. 1. — Structures of (a) temperature and (b) density of the gas within the He core of the 
SN 2006jc model (red) at day 50 (solid lines) and day 200 (dashed lines) after the explosion. 
For comparison, shown are those for the SN II model (blue) with M pr = 25 M and E 5 i = 1 
(Umeda & Nomoto 2002), which has a He core of 5.5 M comparable to that of SN 2006jc. 
The mass coordinate is normalized by the mass of the He core. 
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Fig. 2. — (a) Condensation times and (b) average radii of dust grains formed in the ejecta 
of SN 2006jc as a function of mass coordinate. 



-22 - 



2000 



200day 



C (60days) 



^ 1500 

s 

b" 

2 iooo 

cd 

CD 

6 

CD 



FeS 



C (200days) 



500 



MgSiCL Mg 2 Si0 4 



SiO ? 



A1 2°3 



^LmsO 



_i i i_ 



_i i i_ 



3 4 5 6 

mass coordinate; M r (M Q ) 



Fig. 3. — Temperature of each dust species formed in the ejecta of SN 2006jc at 200 days 
after the explosion. The dashed line depicts the temperature of C grains at 60 days. 
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Fig. 4. — Spectral energy distribution by thermal radiation from the newly formed dust. 
The magenta line is the spectrum obtained by adopting the mass of dust Mij from the dust 
formation calculation. The red line is the best fitted spectrum obtained by taking the dust 
masses and upper mass limits M2J given as in Table 1, for which the contribution from each 
dust species is depicted by the dashed lines. The cyan symbols are the photometric data at 
~200 days by AKARI (Sakon et al. 2007) and MAGNUM (Minezaki et al. 2007). 
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Fig. 5. — Spectral energy distribution derived in the presence of CS carbon grains (thick 
red line), where the mass of newly formed and pre-existing C dust is 2.24 x 10 -4 M & (black 
dashed line) and 1.4 x 10 -3 M Q (black dotted line), respectively. The mass of other newly 
formed dust is the same as that in Figure 4. The cyan symbols are the photometric data at 
~200 days by AKARI (Sakon et al. 2007) and MAGNUM (Minezaki et al. 2007). 
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Table 1. Mass of Each Dust Species and References for Optical Constants 



dust species 


Mij (M Q ) 


M 2d (M ) 


references 


C 


0.701 


5.6 x 10~ 4 


1 


A1 2 3 


0.008 


<0.008 


2 


MgSiOs 


0.157 


<0.157 


3 


Mg 2 Si0 4 


0.082 


<0.082 


4 


Si0 2 


0.229 


<0.229 


5 


MgO 


0.010 


<0.010 


6 


FeS 


0.067 


0.002 


4 


Si 


0.196 




7 


Total 


1.450 


<0.489 





Note. - The mass of grain species j obtained by 
the dust formation calculation is denoted by Mij. The 
mass or upper mass limit of the dust necessary for re- 
producing the SED observed at 200 days is denoted by 
M 2 j, where Si grains are excluded because their con- 
densation time could be much later than 200 days (see 
text). References for the optical constants used in the 
calculation are shown in the last column. The optical 
constants of MgSi0 3 at A < 0.3 /xm are replaced by that 
of Mg 2 Si0 4 (see Hirashita et al. 2008). 

References. — (1) Edo (1983); (2) Toon et al. (1976); 
(3) Dorschner et al. (1995); (4) Semenov et al. (2003); 
(5) Philipp (1985); (6) Roessler & Huffman (1991); (7) 
Piller (1985) 



